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Within A Multi-Phase Transport model, we investigate decorrelation of event planes over pseudo- 
rapidity and its effect on azimuthal anisotropy measurements in relativistic heavy-ion collisions. 
The decorrelation increases with increasing n gap between particles used to reconstruct the event 
planes. The third harmonic event planes are found even anticorrelated between forward and back- 
ward rapidities, the source of which may root in the opposite orientation of the collision geometry 
triangularities. The decorrelation may call into question the anisotropic flow measurements with 
pseudorapidity gap designed to reduce nonflow contributions, hence the hydrodynamic properties 
of the quark-gluon plasma extracted from those measurements. 
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Introduction Relativistic heavy-ion collision data 
indicate that a strongly interacting quark-gluon 
plasma(QGP) is formed where the relevant degrees of 
freedom arc quarks and gluons |l|. In a non-hcadon 
heavy-ion collision, the geometrical overlap region- 
wherc interactions take place between the participant 
nucleons-is elliptic on the transverse plane perpendicular 
to the collision axis. Due to interactions, the high energy 
density and pressure built up in the center of the collision 
region power an anisotropic expansion and collective 
motion of the QGP. This results in an elliptical distri- 
bution in the final-state particle azimuthal distribution, 
called elliptic flow Q. The measured elliptic flow is so 
large that hydrodynamical descriptions are applicable 
and the shear viscosity to entropy density ratio (r//s) 
cannot be much larger than the conjectured quantum 
low limit of 1 /47T Q| ■ Similar phenomenon has been also 
observed in a gas of cold Fcrmionic lithium-6 atoms, a 
system very different from the QGP, where a magnetic 
field is used to induce strong, resonant interactions 

It was not realized until recently || that there can be 
a triangular shape component in the transverse overlap 
region in the configuration space because of fluctuations 
in the nucleon distributions inside nuclei. This triangu- 
larity can result in a third harmonic in the azimuthal 
distribution of final particles, called triangular flow. Hy- 
drodynamical calculations indicate that triangular flow 
is more sensitive to the rj/s [f|. Although the initial con- 
figuration space information of the overlap region is not 
directly observable, their footprint is contained in the 
final-state particle correlations. This is analogous to the 
nonuniform cosmological microwave background as a re- 
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suit of the primordial density fluctuation of the universe 
by gravitational interactions Q • 

Heavy-ion experiments measure anisotropic flow via 
final-state particle correlations. For example, one con- 
structs an event plane (EP) to be the maximum par- 
ticle emission direction, as a proxy for the participant 
plane-the minor symmetry axis of the nuclear overlap 
region in the configuration space. One then correlates a 
test particle with the event plane to measure anisotropic 
flow [1]. As such, the measured anisotropy is contam- 
inated by other particle correlations, generally referred 
to as nonflow Many nonflow correlations are short 
ranged, such as resonance decays and jet-corrclations. 
Thus, to reduce nonflow contributions, one often applies 
a pseudo-rapidity (77) gap between the particles used for 
EP construction and those used for anisotropy measure- 
ments. The basic assumption is that the participant 
plane is the same for all pscudorapiditics. This may not 
be true as first pointed by Petersen et ai. fiof . This is 
because, unlike the reaction plane which is unique in a 
given event, participant planes can be different in differ- 
ent phase-space regions of the same event. In this study, 
wc investigate the event plane correlations in pseudora- 
pidity in a theoretical model and find they are indeed 
different. In other words, it may not be reliable to mea- 
sure anisotropic flow using the nonflow-reducing 77-gap 
method. This finding, if also true in real data, would 
have important implications in terms of the QGP prop- 
erties one extracts by comparing data to hydrodynamic 
calculations. 

Analysis Method We use the AMPT (A Multi-Phase 
Transport) model with string melting for our study. 
There are four main components in AMPT: the ini- 
tial conditions, parton interactions, hadronization, and 
hadron interactions. The initial conditions arc obtained 
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from the HIJING model [ll|, which includes the spatial 
and momentum information of minijet partons from hard 
processes and strings from soft processes. The time evo- 
lution of partons is then treated according to the ZPC 
parton cascade model [12]. After parton interactions 
cease, a combined coalescence and string fragmentation 
model is used for the hadronization of partons. The 
scattering among the resulting hadrons is described by 
a relativistic transport (ART) model [Hj] which includes 
baryon-baryon, baryon-meson and meson-meson elastic 
and inelastic scatterings. 

The azimuthal anisotropy is usually characterized by 
the Fourier coefficients M: 
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where <f> is the particle azimuthal angle and ty n is the 
n-th harmonic event plane angle. Note "I^ is not nec- 
essarily the reaction plane (the plane defined by beam 
direction and impact parameter) due to event-by-event 
fluctuations. 

In AMPT model ^> n can be calculated in coordinate 
space by [f| 

1 7T 

* r = -atan2((r 2 $m.(n<j> part )) , (r 2 cos(n</> part ))) + -, 
n n 

(2) 

where r and 4> pa rt are the polar coordinate positions of 
each parton. However, the coordinate space information 
is not accessible by experiment. The event plane is in- 
stead constructed from measured particle momenta by 

*p = -atan2((sin(n0)), (cos(n0)», (3) 
n 

where <j> is the azimuthal angle of the particle momen- 
tum. We therefore also study the momentum space "J^ 
in AMPT. 

Due to the finite multiplicity of constituents, the con- 
structed event plane is smeared from the true one by 
a resolution factor. The observed anisotropy parameter 
needs to be corrected by the corresponding event-plane 
resolution as 
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The resolution factor can be obtained by the sub-event 
method with an iterative procedure Q. Because of the 
large initial parton multiplicity, our calculated 0^ and 
^3 are nearly unity, even in the most forward or back- 
ward r\ range (3.5 < \r\\ < 4). However, the resolution 
on the final-state momentum space event planes deviate 
significantly from 1. For 20-50% centrality Au+Au colli- 
sions, Ml (01) decreases from 0.302 (0.053) at 77=0-0.5 
to 0.072 (0.008) at r?=3.5-4. 

Results and Discussions Figure Ha) shows the 
event probability distribution in the difference of the 
configuration space event-plane angles and ^3) con- 
tracted from initial parton transverse coordinates in the 
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FIG. 1: (Color online) Event probability distribution in the 
difference of event planes constructed by (a) initial parton 
coordinates and (b) final particle momenta in the forward and 
backward r\ ranges for 20-50% centrality in Au+Au collisions 
at V /5^F = 200 GeV from the AMPT model. 
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FIG. 2: (Color online) Pseudo-rapidity gap dependence of 
the correlation strength (cos[?i(*I' n {77} — ^ n{—rj})]) corrected 
by the corresponding event plane resolutions for 20-50% cen- 
trality in Au+Au collisions at ^/sjvjv = 200 GeV from the 
AMPT model. 

forward and backward rj ranges. Sharp peaks are ob- 
served when are constructed not far from midrapidity 
for both the elliptic and triangular harmonic planes. This 
indicates a strong correlation between ^{^j and ^{—77}. 
With increased r\ gap between the forward and backward 
regions, the probability distribution broadens for in_ 
terestingly, the 'Fg from forward and background 77 with 
large r\ gap are anticorrelated on average. 

The probability distribution in event-plane angle dif- 
ference constructed from final-state charged particle mo- 
menta is shown in Fig. [Tf b) . Similar to the results in 
panel (a), the event plane angle correlations are relatively 
stronger close to midrapidity and become weaker when 
the r\ gap increases. Especially, the triangular harmonic 
plane at very forward and backward 77 are random with 
respect to each other. The event-plane correlations from 
final-state particle momenta is weaker than those from 
the initial-state parton configurations. This is due to the 
worse resolution of the momentum space event plane, a 
direct result of much fewer final-state particles than the 
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FIG. 3: (Color online) Event probability distribution in (a) ^2 and (b),(c) it 3 constructed by initial parton coordinates and 
final particle momenta in different rj ranges for 20-50% centrality in Au+Au collisions at y / sjvjv = 200 GeV from the AMPT 
model. The panel (d) is a cartoon of interacting nucleons from the two nuclei. 



initial partons. 

The decorrclation between the forward and backward 
event planes are, in part, due to finite event-plane res- 
olutions. We therefore show in Fig. [5] the correlation 
strength (cos[n(^ n {rj} — ^ n {— r)})]) divided by the corre- 
sponding event plane resolutions S£ n \r\\ x 3i n {—r]}. The 
correlation is plotted as a function of the 77 gap (A77, de- 
fined as the difference between the centers of the 77 regions 
used for event-plane calculation) ; Ar/ = 2ry for our choice 
of the symmetric 77 ranges. Ideally, when A77 is small 
close to zero, "^{r]} should be the same as ^{—77}. In- 
deed Fig. [2] shows the event- plane correlation approaches 
unity at small A77. However, we find that the correlation 
decreases with increasing A77, falling siginificantly below 
unity. This indicates that the event-plane decorrelation 
is not simply due to the degrading event-plane resolution, 
but physics-the event planes from forward and backward 
pscudorapiditics arc indeed different. The decorrclation 
of ^2 an( i ^2 are similar as function of At;. The decorre- 
lations in the triangular event planes are much stronger 
than those in the elliptic ones. They also have differ- 
ent A77 dependences; The decorrelation magnitudes dif- 
fer between ^3 and "Fg. With large enough ?7-gap, the 
triangular event planes from configuration space arc anti- 
corrclated. 

To illustrate the physics further, we show in Fig. [3ta) 
the event probability distribution in ^2, constructed 
from both the initial parton coordinates and the final 
particle momenta in different 77 ranges, relative to the 
reaction plane angle (fixed at zero). The forward and 
backward vf^ distributions are consistent with each other. 
The ^2 distribution is wider than the A^ distribution in 
Fig. [TJ even though there is decorrelation beyond event- 
plane resolution in A^- This suggests that the recon- 
structed event-plane does not reflect the reaction plane 
but the participant plane; There exists an additional 
fluctuation in the participant plane about the reaction 
plane. In other words, the two non-identical participant 
planes (one at forward and the other at backward pseudo- 
rapidity) are closely correlated, but both deviate signifi- 
cantly from the reaction plane. 

The ^3 distributions are displayed in Fig. [3Jb) . Unlike 



^2; the forward and background ^3 are anticorrelated, 
especially between those with a large A77 gap. As the 
system expands, the property of this anti-correlation in 
the initial geometry is transfered to the final state parti- 
cles in the momentum space, although at a much reduced 
magnitude. This is shown in Fig. [3{c). 

The anti-correlation in ^3 may root in the collision 
geometry. This is illustrated in the cartoon of Fig. [3Jd). 
The nucleons in the overlap region of the forward-going 
projectile nucleus (indicated by the solid sphere) see vari- 
ous thicknesses of the backward-going target nucleus (in- 
dicated by the dashed sphere). The projectile nucleons 
indicated by the dark dots suffer fewer collisions than 
those indicated by the lighter ones, thus have a relatively 
larger probability to end up in the forward direction. In 
addition there are more of the "dark" nucleons than the 
light ones in the projectile nucleus. Therefore, at forward 
rapidity, the distribution of the struck nucleons is more 
likely to have a upward triangular shape (indicated by the 
solid triangle). Conversely, that at backward rapidity is 
more likely shaped by the oppositely oriented dashed tri- 
angle. Thus, triangularities of the forward and backward 
interaction zones are more likely anticorrelated. 

Experimentally, to reduce nonflow, one often measures 
particle anisotropy by applying a 77-gap between the par- 
ticles of interest and those used for event plane recon- 
struction. The event-plane decorrelation we have shown 
calls this 77-gap method into question. This was first 
pointed out by Petersen et al. [10| . Here we study the ef- 
fect of event-plane decorrclation on anisotropy measure- 
ments. In Fig. 21 we investigate vi and U3 of charged 
particles within < 0.5 with respect to \E ,r and *f! p 
from different 77 ranges as functions of the 77 gap. The 
v n is already corrected by the corresponding event-plane 
resolutions. As the r\ gap increases, the magnitudes of 
V2 and V3 decrease. The decrease is, in part, due to the 
decorrelation between the event plane at forward 77 and 
the event plane at |t;| < 0.5. The decorrelation increases 
with increasing A77. This study suggests that one may be 
required to measure particle anisotropy using the event 
plane reconstructed from the same 77 region as the parti- 
cles of interest. 
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FIG. 4: (Color online) Pseudo-rapidity gap dependence of (a) 
V2 and (b) V3 with respect to 9 r and from different 77 
ranges with 3 mb parton cross section for 20-50% centrality 
in Au+Au collisions at ^/s^F = 200 GeV from the AMPT 
model. 

The decreasing behavior of v n with A77 is found for 
both ^> r n and However, the interpretation is compli- 
cated for v n using due to nonflow. Nonflow is gen- 
erally also a decreasing function of A77. The decreasing 
trends of v n with respect to is therefore a result of 
combined effects from the reduced non-flow and decorre- 
lation in the event plane angles over 77 gap. 

Is the difference in v n from ^> r n and from "F^ en- 
tirely due to nonflow? The answer is no, because the 
event plane decorrelations arc different for ^3 and from 
^3, which must result in a correlation (cos[3( v I'3{r?} — 
^3{?7})]) varying with 77. Even with the similar decorre- 
lations for ^2 an d ^2 as a function of 77, we found that 
(cos[2(<3>2{?/} — ^i 7 ?})]) varies with 77. In other words, 
the difference in v n between ^> n and "F^ is a combined 
effect of nonflow and (cos[n(^ n {ri} - $£{77})]). The fact 
that ^{77} and ^{77}, after resolution corrections, are 
not the same may suggest that the final state particles 
in a given 77 region are not solely determined by the ini- 
tial partons from the same 77 region. This is in fact not 
completely unexpected but the discussion of its physics 
is beyond the scope of this paper. 

Conclusions By utilizing AMPT model with string 



melting, we have studied the correlations between event 
plane azimuthal angles reconstructed from particles in 
the forward and backward pcudo-rapidity regions in 
Au+Au collisions at ^/smn = 200GeV. The event planes 
recontructed from the initial parton coordinate space and 
the final particle momentum space are both studied. The 
event plane correlations are corrected by the event-plane 
resolutions. The resolution corrected event-plane corre- 
lation is found to weaken as the forward-backward 77-gap 
increases. This indicates that the event planes from dif- 
ferent 77 ranges differ. The decorrelation in the triangu- 
lar event plane is significantly stronger than that in the 
elliptic one. Particularly, the initial parton coordinate 
^3 exhibits an anti-correlation between very forward and 
backward 77. The anticorrelation may root in the collision 
geometry-the triangularity of the forward going partici- 
pants is opposite to that of the backward going ones. 

The decorrelation may call into question the 
anisotropic flow measurements with pseudorapidity gap 
designed to reduce nonflow contributions. Because of 
the event-plane decorrelation, azimuthal anisotropy mea- 
sured by the event planes from different 77 ranges yield 
different results. By applying a large 77 gap between the 
anisotropy measurements and the event-plane recontruc- 
tion, one can reduce the nonflow contributions, but at the 
same time increase the event-plane decorrelation effect. 
It may therefore not be a reliable strategy to apply a large 
77 gap to measure anisotropic flow. The large 77-gap may 
result in an under-measured anisotropic flow because of 
the event-plane decorrelation. If this were true in real 
data, it could have important implications to the QGP 
properties extracted by comparing data to hydrodynamic 
calculations. In particular, one may speculate that the 
extracted shear viscosity to entropy density ratio could 
be overestimated due to an under-measured flow, which 
in turn could suggest that the QGP may be more perfect 
as we thought it was. 
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